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ABSTRACT 

Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star¬ 
forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With 
the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such 
an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure 
the molecular gas content of seven high-redshift (z ^ 1.6) starburst galaxies with the Atacama Large 
(sub-)Millimeter Array and IRAM/Plateau de Bure Interferometer. Our targets are selected from the 
sample of Herschel far-infrared detected galaxies having star formation rates (^300-800 Mq yr“^) 
elevated 4x) above the star-forming ‘main sequence’, and included in the FMOS-COSMOS near- 
infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. We detect GO emission 
in all cases at high levels of significance, indicative of high gas fractions (~30-50%). Even more 
compelling, we firmly establish with a clean and systematic selection that starbursts, identified as 
main-sequence outliers, at high redshift generally have a lower ratio of GO to total infrared luminosity 
as compared to typical main-sequence star-forming galaxies, although with a smaller offset than 
expected based on past studies of local starbursts. We put forward a hypothesis that there exists a 
continuous increase in star formation efficiency with elevation from the main sequence with galaxy 
mergers as a possible physical driver. Along with a heightened star formation efficiency, our high- 
redshift sample is similar in other respects to local starbursts such as being metal rich and having a 
higher ionization state of the interstellar medium. 

Subject headings: galaxies: ISM — galaxies: high-redshift — galaxies: starburst — galaxies: star 
formation 
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1. INTRODUCTION 


It has recently come to light that the growth of 
galaxies may be less erratic than previously thought. 
The existence of a tight relation between the stel¬ 
lar mass (M*) and star formation rate (SFR) of star- 


Noeske et al.ll2007l: IDaddi et al.ll200' 

?i: lElbaz et all 120071 

Whitaker et al.ll2012l: IKashino et al.l 

2013HSDeagle et al.l 

2014D is indicative of a quasi-equilibrium being main- 


tained betwee n gas supply, SFR, and gas expulsio n from 
galaxies ('e.g.. lBouche et al1l201Cll : lLillv et al.ll2013ll . This 
SFR-M* relation evolves strongly with look-back time, 
paralleling the global evol ution of the cosmic SFR den sity 
from the present to z ~ 2 (|Madau fc Dickinsonll20l3) . It 
is becoming clear that such a decline in SFR is in re¬ 
sponse to diminishing gas reservoirs (e.g., iTacconi et al.l 
l20I3t iScoville et al.l 12014 iSantini et al.l 120141) . The effi- 
ciency of star formation {SFE = SFR/Mgas), i-e., the 
efficiency at which gas is being converted to stars, is 
remarkably similar across a wide range of cosmic time 
within the global star-forming population, usually with 
gas mass Mgas in ferred from the CO molecular line lumi¬ 
nosity L'co ITacconi et al.ll2oM 120131 : iDaddi et al.l 
I2010ijl . 


However, outliers from the MS are known to exist, with 
very high specific SFR (sSFR) compared to normal MS 
galaxies, such as local ultra-lumin ous infrared galaxies 
(ULIRG. ISanders fc Mirabel 19960 . At z ^ 2, these MS 
outliers represent ^ 2% of the star-forming population 
and contribute only a moderate fraction 10%) of th e 
global co-moving SFR-density (iRodighiero et al.l 120111) . 
Yet, they are still likely to play an important role in 
galaxy evolution, as a high fraction of galaxies may ex¬ 
perience such a star-bursting event. Therefore, it is im¬ 
portant to understand the physical conditions that can 
lead to such a surge in SFR. 

Local starbursts have a higher SFE compared to MS 
galaxies. Estimates of their molecular gas content from 
their CO luminosity indicate tha t a given gas mass i s 
capable of producing higher SFRs (|Solomon et aklUOQ^ . 
and a bi-modal distribution of SFEs has b een proposed 
for the general MS and o utlier population (iDaddi et al.l 
l2010bl : [Genzel et 311120101) . While there are uncertainties 
with respect to the appropriate conversion of to H 2 
gas mass (aco), an offset with respect to MS galaxies 
is already apparent when considering the observed quan¬ 
tities Lqq and total infrared luminosity (Ltir)- The 
physical mechanism responsible for this enhanced SEE 
is not entirely clear and one cannot exclude that sam¬ 
ple selection may be responsible for the lack of galax¬ 
ies with intermediate SFEs. While CO detections have 
been obtain ed for high-redshift submm-selected gal axies 
(SMGs; e.g.. lGreve et al.l[2005HBothwell et al.ll2oi3[) an d 
star-forming radio galaxies fSERG: iCasev et all 120111) . 
these samples do not cleanly distinguish starburst out¬ 
liers from the MS population and submm selected sam¬ 
ples could be biased towards gala xies with lower dust 
temperatures (|Magnelli et aLllMTll . Therefore, it is im¬ 
perative to study the CO emission of starbursts securely 
elevated above the MS at their respective stellar mass 
and near the peak epo ch of star formation (z 2). 
IRodighiero et akl (120 Ilf) achieve such a clear distinction 
between MS and MS-outlier galaxies using large photo¬ 


metric samples with multi-wavelength coverage, includ¬ 
ing the bulk of the EIR emission. Our targets then result 
in a homogenous sample that does not have contaminants 
such as those inherent with other selections. 

To address this issue, we report on interferometric ob¬ 
servations of the molecular transition ^^CO 2-1 and ^^CO 
3-2 for a sample of seven MS outliers using both the 
Atacama Large Millimeter Array (ALMA) and IRAM 
Plateau de Bure Interferometer. Our sample is extracted 
from the 246 MS outliers (with sSFRs > 4x above the 
mean of the MS; Figure [T]) at 1.4 < z < 2.5 identi- 
fied with Herschel obse rvations over the COSMOS field 
(|Rodighiero et al.ll2011h that effectively cover the peak 
FIR emission. Among these Herschel-detected galax¬ 
ies, those with photometric redshifts between 1.4 and 
1.7 are observed throu gh our Subaru Intensive Project 
(jSilverman et al.1 [20151 1 with the FMOS near-IR multi¬ 
object spectrograph. These spectra provide a detection 
of key diagnostic emission lines (i.e., Ha, H/3, [NHJA6585, 
[OHIJA5008) used to measure accurate spectroscopic red- 
shifts, SFRs, and metallicities. Six targets, labeled in 
Figure [TJ are extracted from this sample with an ad¬ 
ditional galaxy at a slightly higher redshift (PACS-282; 
z = 2.19) wher e the redshift com es from the zCOSMOS- 
Deep program ((Lillv et al.ll2007ll . 

Our confidence in these galaxies being ‘outliers’ in 
the SFR—M* plane is based on the exquisite multi¬ 
wavelength coverage of the COSMOS field. Stellar 
masses are derived from SED-fitting using Hyperz with 
stelly population synthesis models (jBruzual fc Charlotl 
1200311 at the respective spectroscopic redshift. The de¬ 
blending of detections in Herschel (or 70/rm Spitzer) im¬ 
ages relies upon Spitzer MIPS 24^m priors. SFRs are 
derived from Lttr. an integral of the best-fit model SED 
(|Draine fc Lill200'^ from 8 - 1000 /rm using the Spitzer 
(24 p,m), Herschel PACS (100 and 160 /rm) and SPIRE 
(250, 350 and 500 /im) band s, thus accouri t ing fo r the ob¬ 
scured star formation as in iMagdis et al.l (1201211 . More¬ 
over, four of the seve n galaxies have 1.4 GH z radio de¬ 
tections at > 5cr level (iSchinnerer et al.ll20in[l and radio- 
based SFRs are consistent with being above the MS. The 
presence of strong AGN within our sample is ruled out 
by the lack of individual X-ray detections by Chandra. 
Even if an obscured AGN were present, the contribu¬ 
tion to LTfB wou ld be negligible (|Rodighiero et al.ll2011l : 
iPozzi et al.ll2012 1. Throughout this work, we assume 
Flo = 70 km s“^ Mpc“^, Ha = 0.7, Hm = 0.3, and 
use a Salpeter IMF for SFRs and stellar masses. 

2. OBSERVATIONS OF CO EMISSION AT HIGH-Z 

ALMA observations of five galaxies were carried out 
in Cycle 1 (Project 2012.1.00952.S) using 25-30 12m an¬ 
tennas and the band 3 receiver tuned to intermediate 
frequencies between 86.395 - 94.888 GHz with a reso¬ 
lution of 488 kHz and spectral bandpass of 1.875 GHz. 
Two ‘scheduling blocks’ (SBl and SB2) were defined to 
accommodate the requirements to detect CO 2-1 within 
a single sideband for each of our targets (SBl: PACS- 
819, 830; SB2: PACS-325, 299 and 867). On-source ex¬ 
posure times were 32 minutes (819 and 830), 2.3 hours 
(299 and 325) and 3.8 hours (867) that met our request 
based on predictions of CO emission from their IR lumi¬ 
nosity and an assumed ratio Lco/Ltir (3x lower than 
MS galaxies). The antenna configurations had baselines 
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Table 1 

Sample data and CO measurements^ 



ID 

RA 

(CO) 

Dec 

(CO) 

zco'" 

? c 

-^spec 

log mN 

M0 

SFR® 

Jy km s~^ 

f' g 
^CO 

Av*' 

CO 

size (") 

Ghs 

fraction 

299 


09:59:41.31 

02:14:42.91 

1.6483 

1.6467 

10.44 

554±40 

0.67±0.08 

10.44 

590 

<2.4 

0.52 

325 


10:00:05.47 

02:19:42.61 


1.6557 

10.29 

139±10 

0.28±0.06 

10.07 

764 


0.32 

819 


09:59:55.55 

02:15:11.70 

1.4451 

1.4449 

10.61 

783±18 

1.10±0.07 

10.55 

592 

0.34±0.08 

0.34 

830 


10:00:08.75 

02:19:01.90 

1.4631 

1.4610 

10.86 

517±12 

1.18±0.10 

10.59 

436 

0.97±0.17 

0.46 

867 


09:59:38.12 

02:28:56.56 

1.5656 

1.5673 

10.67 

358±34 

0.46±0.04 

10.24 

472 

<1.5 

0.29 

282 


10:00:01.54 

02:11:24.27 

2.1869 

2.1924 

10.88 

581±42 

0.75±0.12 

10.44 

660 

<3.4 

0.28 

164 


10:01:30.53 

01:54:12.96 

1.6481 

1.6489 

>10.28 

358±8 

0.61±0.11 

10.40 

894 

<4.8 

0.52 


^ The first five targets are observed with ALMA while the remaining two with PdBI. 

^ (Tz — 0.0003-0.0004. An accurate CO centroid was not obtained for 325 thus the Ha redshift was used for the CO flux measurement. 
^ Spectroscopic redshifts are based on Ha with the exception of and have errors <y /{ i -\- z ) — 1-8 X 10~^. 

^ ctm ~ 0.07 dex jllbert et al.|[^015ll 
® Units of Mq yr“^. 

^ All are based on the CO 2-1 transition with the exception of 282 (3-2). 

S Log base 10; units of K km s~^ pc^ 

^ Full width of the CO line at zero intensity in units of velocity (km s~^). 


1.4<z<1.7 



Figure 1. SFR versus M* for galaxies with 1.4 < 2 : < 1.7 in 
COSMOS. Galaxies with CO observations are shown in blue. For 
reference, we plot star-forming galaxies that denote the MS (small 
dots). Herschel-detected galaxies include our sample in blue and 
those shown in red. 

up to ^400m (^200m) for SBl (SB2) thus resulting in a 
beam size (FWHM) of 1.3" x 1.0" (4.5" x 2.0"; 5.4" x 2.7" 
for PACS 325 only). The images of PACS-819 (830) have 
a flux sensitivity level of 0.12 (0.10) mJy beam“^ over 
a bandwidth of 400 km s“^ while the deeper observa¬ 
tions reach 0.073 - 0.14 mJy beam“^ for the same spec¬ 
tral window. Standard targets were used for flux (Mars, 
Ganymede), bandpass, and phase calibration. 

We measure the CO flux by fitting the data in UV 
space with the GILDAS task ‘uvfit’ available in the 
MAPPING package. The fit was performed using a point 
source model or a circular Gaussian for resolved sources 
that returns the centroid, de-convolved source size, and 
integrated flux. We find good agreement between fluxes 
returned from GILDAS with those using ‘imfit’ in the 
image plane with CASA. Additionally, we resolve the 
emission for PACS- 819 and 830; therefore, a measure of 
the source extent is reported. 

To increase the sample size, interferometric measure¬ 


ments were obtained with IRAM/PdBI for PACS-282 
(CO 3-2) and PACS-164 (CO 2-1). For PACS-282 
(PACS-164), an integration time of 6.2 (12.9) hours, us¬ 
ing the 3mm band in configuration D, was achieved that 
resulted in a limiting Icr sensitivity level of 0.14 (0.11) 
mJy beam“^ for a bandwidth of 240 (180) MHz. CO 
fluxes are measured with GILDAS. In Table [TJ we list 
the source properties and CO measurements. 

3. RESULTS 

We detect CO emission in all seven targeted galaxies 
with integrated flux densities (Ico) ranging from 0.28 
to 1.18 Jy km s“^. All but one have a high level of 
significance (> 5a) and those resulting from ALMA ob¬ 
servations are above the 8 ct level with the exception of 
PACS-325 {S/N = 4.7). In Figure^ w e display the opti¬ 
cal H ST/ACS F814W image cutouts (jKoekemoer et al.l 
l2007f) with CO and 3.6/rm emission overlaid as contours. 
Using a large beam for the majority of the sample, the 
CO emission is essentially unresolved with the exception 
of two observations taken at higher resolution (^ 1"). 
PACS-819 has half of the emission coming from a re¬ 
gion of 3.2±0.8 kpc, while PACS-830 is more extended 
(9.Oil.6 kpc). This may suggest that there is diver¬ 
sity in the siz e distribution of molecular gas at high 
redshift (e.g . . Ilvison et Tl 120111 : IRiechers et al.l 120111 : 
Hodge et al. ll2?)lT dissimilar to local starbursts (e.g., 
Scoville et al. Ipli). 

Upon close inspection of the maps in Figure [21 it is 
evident that the centroid of the CO emission is not always 
coincident with the brightest regions of UV emission as 
seen in the HST/ACS image (i.e., PACS-164, PACS-830, 
PACS-867). This is likely evidence that a fair fraction 
of the star formation is obscured as supported by the 
improved alignment of the CO emission with the peak 
infrared emission detected by Spitzer in both the IRAC 
and MIPS channels, and clear associa tion with highly 
signif icant radio emission at 1.4 GHz (ISchinnerer et al.) 

HUg). 

We note that the Herschel far-IR (and possibly the CO) 
emission of PACS-164 is likely the sum of two separate 
components seen in the IRAC image. The peak of the 
CO emission, detected by IRAM, is located somewhat 
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Figure 2. HST/ACS F814W images of our seven starbursts. CO emission is shown as solid (dashed) blue contours in positive (negative) 
steps of 1(T starting at 2 (-2) times the noise level (where arm3=0-i7, 0.12, 0.11, 0.05, 0.12, 0.14 and 0.11 mjy beam“^ starting from left 
to right with 282 and then 819). Spitzer/IRAC 3.6/rm detections are marked with red contours to illustrate the association of peak CO 
emission with their IR counterparts. CO beam size is indicated in the lower left corner. Small black circles show the placement of FMOS 
fibers. For PACS-164 the R-band image is shown to further illustrate the heavily-obscured nature of the CO-emitting region. Angular 
sizes of the cutouts varies. 


between two IR sources, as made particularly clear by 
the itl-band image also shown in Figure O and slightly 
elongated possibly indicating a contribution from both 
sources. The FMOS spectrum that provides the redshift 
{z = 1.650) refers to the western IR component that is co- 
spatial with the UV-bright source seen in the HST image. 
With a redshift of the CO emission {z = 1.647) very 
close to that of the FMOS source, this system appears 
to be in the early stages of a merger given the projected 
separation of 18.7 kpc. The lower limit for the stellar 
mass of the system as reported in Table 1 refers to the 
mass of the object observed with FMOS, and encircled 
in Figure 2. 

3.1. CO-to-IR ratio 

Our primary interest is to determine whether high- 
redshift starbursts convert gas into stars with a higher 
efficiency than MS galaxies at these epochs. To begin 
with, we use the observed CO to IR luminosity ratio 
Lqo/^tir rather than derived quantities, e.g. the gas 
mass, that depends on the CO-to-H 2 conversion factor 
(«Co); or gas mass surface density, which would fur¬ 
ther require CO-size information for which we only have 
for 2/7 galaxies in our sample. However, we note that 
Lqo/-^tir oc Mgas/SFR (modulo aco), hence this lu¬ 
minosity ratio is a fair proxy for the gas depletion time 
which in turn is the inverse of the SFE. 

In Figure (HJi, we plot ^CO as a function of Ltir, in¬ 
dicative of the obscured SFR, and include both low- and 
high-redshift galaxies with measurements availa ble in the 
literature and compiled in iSargent et al.l (j2014[) . All line 
luminosities are converted to CO (1-0) using values of 
0.85 and 0.7, respectively for CO (2-1) and CO (3-2) 
(Daddi et al. 2015). Different excitation corrections 
than those applied here would not impact our results 
since changes would either be insignificant (5 - 10% if 
CO (2-1) transitions were generally more excited than 


assumed here) or affect only one of our sources (PACS 
282, for which apply ing the CO(3-2)/CO(l-0) ratio of 
iBothwell et al.l (I20R^ would result in a 35% difference). 

All galaxies in our sample have Lqq/Ltir below the 
well-defined correlation for MS galaxies. These ob¬ 
servations indicate the existence of an offset in the 
Ln ri/ATTR ratio f o r high -redshift starbursts as found 
bv iSolonion et al.1 (|1997t) for local ULIRGs. In Fig¬ 
ure Et, we plot these quantities, normalized to the mean 
value of galaxies on the MS. In addition, on the ab¬ 
scissa we replace Ltir with its implied specific SFR 
(sSFR=SFR/M*), also normalized to the MS value. We 
measure a median value of j ^ L'qo >ms= 0.60 ± 
0.04 for our sample. This represents an offset (1.7±0.1x) 
from the MS that is smaller than in local samples of 
starbursts (~ 3x; dashed line in Fig. [SJa). Along these 
lines, the starbursts in our sample have higher L'qq than 
expected by ^ 0.2 dex (given their excess sSFR rela- 
tive to MS g a laxies ) based on the empirical model of 
ISargent et al.l (j2014D . shown in Figure!^. 

These results may be indicative of a more continu¬ 
ous range in SFE, following the trend for MS galaxies 
(Fig. [3l^) and opposed to the notion of a second mode 
of star-formation operating at higher efSciency distinct 
from MS galaxies (iDaddi et al.ll2010bl : lGenzel et al.ll2010l 
I2015D . However, any statement on differences in SFE re¬ 
quires a factor aco as further addressed below. Even 
with continuity in these parameters, there could still ex¬ 
ist an underlying bi-modal physical framework for star 
formation as demonstrated with e mpirical models shown 
in Eigure[3l) and fully explained in ISargent et al.l (j2014D . 

3.2. Gas masses, aco, cind SFE 

With PAGS-819 and PAGS-830 having marginally re¬ 
solved GO emission, we can estimate t he total gas mas s 
using the dynamical mass method (e.g.. lTan et alllMl . 
The dynamical mass (M<jyn) within the half-light radius 
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SFR [M^/yr] 

Figure 3. (a) CO luminosity versus total infrared luminosity. (6) CO luminosity and specific star forma tion rate with each n ormalized 

to the typical values of MS galaxies. Empirical model curves and Icr errors (grey region) are described in ISargent et all d2014l') . (c) Gas 
depletion time (rdei = SFE~^; units Gyr) versus SFR. Colored symbols with error bars show our starburst sample with two (819 and 
830) in bright red h aving dynamical mass measurements. Grey symbols represent published samples with chqo estimates as compiled in 
ISargent et al.l II2014I1 and shown in all panels. 


(ri/ 2 ) is given for the spherically symmetric case by: 

5a^ri /o 

Mdy„(r < ri/ 2 ) ~ —-— (1) 

where (t= Ar;co/2-35 is the velocity dispersion, Avco is 
the line width - FWHM, and G is the gravitational con¬ 
stant. The gas mass is then derived by subtracting from 
the dynamical mass half of the stellar mas s and a dark 
matte r component (Mdm = 0.25x Mdyn', iDaddi et al.l 
1201 nail as follows. 

-Viyn = 0.5 X (ikf* -|- Mgas) + < f'l/2)- (2) 

With this method, the gas mass for PACS-819 (PACS- 
830) is estimated to be 2.0 ± 1.2 (6.2 ± 2.2) x 10^° Mq, 
based on their half-light radii, respectively 1.4 ± 0.3 kpc 
(4.1 ± 0.7 kpc) and velocity dispersion a, respectively 
168.0 ± 20.2 km s“^ (148.7 ± 20.8 km s“^), that cor¬ 
responds to a gas fraction /gas = Mgas/(Mgas + M*) 
of 0.34 (0.46). Based on these values, we estimate 


aco = 0.6±0.3 for PACS-819 and 1.6±0.6 for PACS-830, 
both not far from those found for the high-re dshift sub¬ 
millim eter galaxies in the proto-cluster GN20 (jTan et al.l 
l2014f ) and consiste nt with lower conversio n factors for 
galaxies in mergers (|Naravanan et al.ll2011[l as comp ared 
to isolated disk galaxies fsee also lCenzel et al.ll2015ll . 

We roughly estimate the H 2 gas masses from the CO 
line luminosity for the galaxies with unresolved CO emis¬ 
sion using the average (aco = 1-1 Mq / K km s“^ 
pc“^) value of the two estimates of aco given above 
which is in very close agreement t o that found for lo¬ 
cal ii uclear starbursts (aco = 0.8; iDownes fc Solomc^ 
Il998f1 . Based on this conversion factor and keeping in 
mind the uncertainty, the gas masses are in the range 
(1.3 — 6.2) X 10^° Mq. This results in gas frac¬ 
tions /gas = Mgas/(Mgas + M,) between 28 and 52% 
(see Table [T] for individual estimat es), similar to those 
of MS galaxies at the se redshifts (|Tacconi et al.l 120131 : 
iBethermin et al.l l2015f l and in line with cosm ological 
galaxy formation models (|Naravanan et al]l2012h . 
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Figure 4. Rest-frame emission line ratios from FMOS-COSMOS: 
(a) Mass - metallicity relation of our high-z starburst sample (red 
points). For compariso n, we plot the M -Z relation from SDSS 
and FM OS at 2 : ^ 1.6 llZahid et al.ll2014ll . (b) BPT diagram at 
2 : ~ 1.6 dKartaltepe et al.ll2015ll with our CO sample as indicated. 
Individual measurements , shown in grey, repre sent galaxies from 
our larger sBzK sample IlSilverman et af1l2Q15l) . Contours denote 
the region spanned by low-redshift galaxies in SDSS. Theoretical 
curves separating star-forming galaxies and AGN (that evolve with 
redshift) are given. All data points have dzlo- errors. 


Based on these results, an increase in the SFE is likely 
responsible for their elevation above the MS. In Figure[21:, 
we illustrate this scenario by plotting the gas depletion 
timescale {jdei = SFE~^; units Gyr) as a function of 
SFR. As described above, the lower CO luminosities at a 
given LtiRj compared to MS galaxies, are indicative of a 
lower gas mass hence shorter gas depletion timescales 
and higher SFEs, in agreement with previous studies 
(|Genzel et al.ll2010l : iDaddi et al.ll2010bll . Despite the re¬ 
maining uncertainties on the assumptions used to derive 
H 2 masses, we see evidence that starbursts in our sample 
have a smaller contrast (i.e., less extreme) in SFE, rela¬ 
tive to MS galaxies, as compared to local ULIRGs and 
the strongest starbursts at high-z. 

As a final check, we can assess whether the values of 
aco derived above are appropriate for our starbursts 
by using the metallicity of the ISM as a proxy, given 
that it is well estab lished that ana is anti-correlated 
with metallicity (e.g ..lArim oto. Sofue fc Tsui imotoll 19961 
iSchruba et al.lDOla iGenzel et al.ll2012lb In Figure HJi, 
we show the ma ss-metallicity relation at z ~ 1.6 

(jZahid et al.l I20l4l based on [NII]A6585/Ha, with our 


sample shown in red. We find that 4 out of 6 galaxies 
(having errors < ±0.2 on the ratio) have high metallici- 
ties, similar to those of local massive galaxies. Therefore, 
our derived values of aco are consistent with local galax¬ 
ies of a similar mass and metallicity, although with the 
caveat that we are not certain whether the metallicity of 
the line-emitting gas is representative of the molecular 
gas producing the CO emission. 

We also examine whether AGN photo-ionization is con¬ 
tributing to the high [Nil]/Ha ratios seen in some of our 
galaxies. Figure® shows the FMOS-COSMOS version 
of the BPT diagram at ^ ~ 1.6. Four galaxies with 
CO measurements have all key diagnostic lines detected 
by FMOS. None of our galaxies fall in the region where 
strong AGN contribution is expected, above the line 
separating AGN fro m star-forming galaxies at z ^ 1.6 
(jKewley et al.l[2?)T^ . While PACS-819 approaches this 
line, it is not detected in X-rays with Chandra. Nev¬ 
ertheless, we cannot rule out the presence of a low-to- 
moderate luminosity AGN, especially because its strong 
[Nil] emission also drives this galaxy above the local MZ 
relation. For reference, in Figure®, we plot the data for 
the larger star-forming galaxy population from FMOS- 
COSMOS (grey data points). It is worth highlighting 
that both PACS-830 and PACS-299 fall within the lo¬ 
cus of the star-forming population, clearly offset towards 
a higher ionization state of the ISM compared to local 
galaxies ([Steidel et al.ll20ll: iKartalteoe et al.ll2ni5h . 

4. CONCLUDING REMARKS 

Our observations of the molecular CO gas content of 
seven galaxies well-above the MS at z ~ 1.6 with ALMA 
in Cycle 1 and IRAM/PdBI establish an offset in the 
^'co — Atir relation for starbursts at high redshift, com¬ 
pared to MS galaxies, although smaller than expected 
from previous studies of low-z starburst outliers. These 
results may be indicative of a continuous distribution in 
SFE at high redshift, as a function of distance from the 
star-forming MS, as opposed to a bi-modal distribution. 

An appealing physical explanation for a decrease in the 
L[,q/Ltir ratio for galaxies well-above the MS is galaxy 
mergers that can lead to rapid gas compression hence 
effectively boosting star formation resulting in shorter 
gas depletion timescales. Several galaxies well-above the 
MS are indeed in a merger phase (iRodighiero et al.ll2011l : 
iWuvts et al.ll201lll . Within the present sample, hints of 
merging come from the presence of multiple UV/optical 
emitting regions (e.g., in PACS-819, PACS-830, PACS- 
867), even being in a kinematically-linked system (PACS- 
164). 
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